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Abstract 

The frequency of a classical periodic system and the energy levels of the corresponding quan- 
tum system can both be obtained using action variables. We demonstrate the construction of 
two forms of the action variable for a one dimensional harmonic oscillator in classical, relativistic 
and quantum regimes. The relativistic effects are considered as perturbative, within the context 
of a non-relativistic quantum formalism. The transition of the relativistic quantum system to 
both classical relativistic and classical non-relativistic regimes is illustrated in a unified framework. 
Formulas for the frequency of a classical relativistic oscillator and the energy eigenvalues of the 
corresponding quantum oscillator for the weak relativistic case are derived. Also studied are the 
non-relativistic and classical limits of these formulas which provide valuable insights on the par- 
allels between relativistic and non-relativistic systems on the one hand and between classical and 
quantum systems on the other. 
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I. INTRODUCTION 



Classical periodic systems can be analyzed elegantly in terms of their action and angle 
variables which constitute a set of canonically conjugate momenta and coordinates. Ac- 
tion variables are proportional to <fpjdqj, where (qj,Pj) are the system's coordinates and 
canonical momenta. For separable systems they are constants of motion in the manner of 
angular momentum and energy. The frequencies of periodic systems can be found using 
the functional relationship between the action variables and total mechanical energy with- 
out requiring a complete solution of the dynamical equations. The French astronomer and 
mathematician Charles-Eugene Delaunay (1816-1872) invented action and angle variables in 
the course of his study of periodicity of lunar motion 1 . These variables assumed importance 
during the early days of quantum mechanics. Lord Rayleigh had shown that in a sinu- 
soidally oscillating system such as a pendulum whose string is shortened slowly, the ratio of 
the energy to frequency, which is directly proportional to the action variable, remains a con- 
stant. Following the language of thermodynamics such motion was referred to as "adiabatic" 
motion. At the first Solvay Conference in 1911, which considered the issues that the new 
quantum ideas introduced into mechanics, it was realized that the adaibatic invariance of 
the action variable in atomic systems, in an environment of slowly varying electromagnetic 
fields, would lead to atomic stability without transitions between states^. Schwarzschild 
introduced into quantum theory the analytical method of employing action variables^. The 
quantization rules of Sommerfeld, Wilson and Ishiwara required that the action variables be 
integer multiples of h to account for the energy spectra of atomic systems.-^. Ehrenfest 
formulated the "adiabatic principle" according to which a slow variation of some parameters 
of a periodic system's Hamiltonian would result in a gradual change in the system's motion 
while maintaining the constancy of the action variable^. With the establishment of the 
wave and matrix forms of quantum mechanics the program of employing action variables 
in the quantum context did not receive significant attention. In the JWKB approximation 
scheme for the determination of bound quantum states the energy eigenvalues are obtained 
by discretizing the action variable. Thus the quantum conditions of the old quantum theory 
can be rigorously deduced as an approximate result in the new quantum theory. 

In 1983 Leacock and Padgett*^ presented a form of quantum mechanics, patterned on 
the classical Hamilton- Jacobi theory and equivalent to the Schrdinger theory, whose focus 
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was a quantum version of the action variable which reduces to the traditional action variable 
in the classical limit. In this formalism the dynamical equation is a postulated quantum 
Hamilton- Jacobi equation for the Hamilton's principal function S which generates a trans- 
formation from the coordinates and momenta (xi,Pi) to angle and action variables (tOj, Jj). 
The bound states of a system are characterized by its quantum action variables assuming 
values which are integral multiples of h in a natural way, and not through arbitrary impo- 
sition as was the case in the Wilson- Sommerfeld scheme. Since the total mechanical energy 
of the system is a function of its action variables it too assumes discrete values for bound 
states. The examination of the dynamical equations of this formalism of quantum theory, 
and the relation between quantum action variables and energy, shows the classical-quantum 
correspondence in a new light. A systematic development of canonical transformations and 
the Hamilton- Jacobi theory in quantum mechanics is found in Ref. 11. Many systems ad- 
mitting bound states have been studied using this quantum Hamilton- Jacobi formalism^^. 
Our objectives here are two fold. We show, using a perturbative approach based on two 
equivalent forms of the action variable, how a weakly relativistic quantum oscillator can 
be treated within the non-relativistic quantum Hamilton- Jacobi formalism to obtain its en- 
ergy eigenvalues that incorporate first order relativistic corrections. Two, we show that the 
dynamics of this weakly relativitic oscillator has the correct non-relativistic classical and 
non-relativistic quantum limits. 

The symbols and notation we use for physical variables have the following meaning. 
The angular frequency of a classical non-relativistic simple harmonic oscillator of mass m 
connected to a massless spring of spring constant k will be ujq = y~- The suffix c refers to 
a classical variable; the absence of this suffix indicates that the quantity being referred to 
is its quantum counterpart. Similarly the suffix r refers to a relativistic variable, and wr is 
the suffix for a variable in the weak relativistic case. 

II. ACTION VARIABLE IN CLASSICAL MECHANICS 
A. Classical Hamilton-Jacobi Theory 

The time evolution of a classical system is governed by its Hamiltonian H which is a 
function of its coordinates x iy the conjugate momenta Pi and the time t. The dynamics of 
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such a system is determined by Hamilton's equations of motion 

dH(xi,pi,t) . dH(xi,Pi,t) 

Xi = ~ , Pi = ~ • (1) 

dpi dxi 

The Hamiltonian of a particle of mass m moving in one dimension under the influence of 
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a potential energy function V(x) is given by H = 2— + V(x). Such a time independent 
Hamiltonian is a constant of the motion and is the total energy E of the system. Thus, 

£ + V(.*) = E. (2) 

Canonical transformations transform one set of coordinate and momentum (x, p) to another 
set (X, P) while preserving the form of Hamilton's equations. One such transformation is 
generated by the function Wc{x,P), whose arguments are the "old" coordinate, x, and the 
"new" momentum, P: 

dW c (x,P) dW c (x,P) 
P= dx ' X= dP ■ (3) 

If this transformation transforms the Hamiltonian into a function only of P, then, using 

©, 

■ dH(P) 



dX 



=>- P(t) = P, a constant, 



dH( P) 

X = dp = ^o, a constant, X(t) = V t + X . (4) 

Thus X and P evolve very simply in time; the former has a linear temporal progress and the 
latter is a constant. The function Wc{x, P), which generates a canonical transformation in 
which the transformed Hamiltonian is independent of the new coordinate X, is the Hamil- 
ton's characteristic function. It is related to the Hamilton's principal function So through 
Sc{x,P,t) = Wc(x,P) — Et, and, for the case of time independent Hamiltonians, satisfies 
the Hamilton- Jacobi equation obtained by using ([3]) in ([2]): 

1 I ' dWAx P)\ 2 

1 OVVc ^^> \ +V (x) = E(P). (5) 



2m \ dx 

The use of this method to solve the dynamical problem involves the following steps: (i) 
Define a suitable new constant momentum P, (ii) Integrate Eq. (J5j) to obtain Wc{x, E(P)), 
(iii) Obtain x(X, P) and p(X, P) using Eq. ([3]), and (iv) Express X and P in terms of the 
initial values x ,p and t. 



B. Action-angle variables and periodic motion 



One particular form of Hamilton- Jacobi theory is especially suited to the study of periodic 
motion. If an inspection of the Hamiltonian indicates that the motion is periodic, then by 
a particular choice of the new momentum P we can evaluate the period of motion without 
obtaining a complete solution of the dynamical problem. The new canonically conjugate 
coordinate and momentum are chosen to be X = w, P = Jq with 

Jc = f Pc{x,E)dx, (6) 



where Pc(%, E), from (T2J), is yj2m[E — V{x)\ and the integral in phase space is performed 
over one cycle of the periodic motion. Jc is the classical action variable and w the angle 
variable. We note that the integral for Jq is the area enclosed in phase space by the path of 
the oscillator's orbit. A new momentum, similar to Jc, will be defined in the corresponding 
quantum formalism and will be referred to as J. Since Jc = Jc(E) we can invert it to 
obtain E = E(Jc)- From Eq. the time evolution of the new coordinate is w(t) = ut + wo 
where the constant "velocity" is 

9H(J C ) dE(J c ) 

It can be shown^ that uj is the angular frequency of this periodic motion. Thus the mathe- 
matical problem of finding the frequency of motion for a periodic system is reduced to that 
of performing the integral ([6]), solving for E to get E(Jc), and evaluating dE/dJc- This is 
a simple and elegant method for evaluating the frequency of a system known to be periodic. 
An equally simple method can be used to obtain the energy eigenvalues of bound states in 
quantum mechanics. 

An equivalent definition of Jc, which is useful for extending the action variable into the 
quantum arena, is 

■Jc = -^j p c (x,E)dx, (8) 

where p c {x, E) is a complex valued function of the complex argument x, and is defined as a 
suitable branch of 

Pc(x,E) = y/2m[E-V(x)]. (9) 
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The turning points x\ and x<i are defined by p c (xi,E) = p c (x2,E) = 0. These are also the 
branch points of p c {x, E) in the complex- x plane. We choose a branch cut connecting X\ and 
X2 along the real axis. p c {x, E) is chosen as that branch of the square root which is positive 
along the bottom of the cut. The counterclockwise rectangular contour C wraps around 
this branch cut. The integral in (jSJ) is performed by deforming the contour C outward 
to the circular contour 7 which lies in an annulus in which p c {x, E) is analytic, expanding 
p c (x, E) in a Laurent series in that annulus and integrating using Cauchy's residue theorem. 
Sommerfeld was the first to employ this contour integral technique in evaluating the action 
variable for the bound states of the electron in hydrogenic atoms. 

An alternate construction of the action variable arises from another canonical transfor- 
mation scheme where the alternate Hamilton's principal function Sc{p,X,t) of the "old" 
momentum and "new" coordinate is used. For conservative systems, the canonical trans- 
formation generated by the alternate Hamilton's characteristic function Wc(p,X), where 
S c (p,X,t) = W c (p,X) - Et, is x = - dWc ^ x \ p = - Wc ^ P x X) . The Hamilton- Jacobi equa- 
tion satisfied by this characteristic function is 

*(-^-)=* ^ 

As in the previous scheme, for periodic systems, we can choose the new coordinate as the 
alternate action variable Jc = — j- § c , %c(p, E)dp, where the clockwise rectangular contour 
C in the complex-p plane surrounds the branch cut along the real axis connecting the 
turning momenta p\ and P2 and the sign of xc{p, E) is chosen positive below this cut. We 
will demonstrate the use of this alternate form of the action variable in both the classical 
and quantum contexts for the weakly relativistic harmonic oscillator. The angle variable, 
while important in describing the state of periodic motion, is not essential for determining 
the frequency (and the energy levels in the quantum context) directly, and therefore will not 
be considered here. 

III. CLASSICAL WEAKLY RELATIVISTIC HARMONIC OSCILLATOR 

The relativistic motion of the harmonic oscillator is governed by the Hamiltonian 
H(x,p) = \Jp 2 c 2 + m 2 c 4 + 7}kx 2 . The total mechanical energy of the relativistic oscilla- 
tor will be referred to as E, and E = E — mc 2 is its mechanical energy in excess of its rest 
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mass energy. The dimensionless energy related parameter we will use is e = We define 
the weak relativistic motion of the oscillator as one characterized by e << 1, and evaluate 
all dynamical variables up to the first order in e. We approximate the Hamiltonian for the 
relativistic oscillator by H = — 8 £ 3c2 + ^f-, retaining only the leading relativistic term in 
kinetic energy, and obtain the oscillator's approximate relativistic frequency using the two 
forms of the action variable discussed earlier. Such a study of the classical oscillator, with 
an approximation of the relativistic Hamiltonian, although unnecessary for frequency deter- 
mination (at least three different series solutions exist for the frequency of a fully relativistic 
oscillator 15 ), is a necessary prelude to the consideration of the corresponding quantum sys- 
tem within the framework of a non-relativistic quantum theory, to be discussed in the next 
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section. Using T for the non-relativistic kinetic energy 2— the correctional relativistic term 
in H can be written as — hT(-^~). The factor is of order e, and can be treated as 
a perturbative parameter in the non-relativistic Hamiltonian. This approximation of the 
Hamiltonian is good for studying motion in the weakly relativistic case. For consistency we 
should retain terms only up to order e 1 in expressions for pcwr, Jcwr and uiwr- 

We first consider the classical action variable in the § pdx form. The classical orbit 
equation is Pc ^ R — P ^ W 3 R 2 = E. Extending this into the complex-x plane the momentum 



is a branch of 

PcW r(x, E) = V2mc^l- ^1 - (e - hx^j . (11) 

Two turning points emerge at x\ and X2, identical in form to those in the fully relativistic 
case, given by —x\ = X2 = y~^§-- These are also two of the branch points of Pcwr{ x , E). 
There are two additional branch points, x^wr and x^wr, arising from the weak relativistic 
correction in the Hamiltonian, given by 




— %3WR — %iWR — \ 1 ~ 7T- (12) 



For e << | these branch points are on the imaginary axis with x^wr farther from the origin 
than x%. We choose one branch cut of Pcwr to connect x\ and x<i along the real axis, and 
two other cuts, each connecting x 3 wr an d x^wr to x = oo along the imaginary axis. We 
choose the sign of pcwr{ x , E) to be positive below the cut joining x\ and x^- Pcwr is 
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analytic for x 2 < \x\ < x^wr and its Laurent series in this annulus can be written as 
Pcwr(x, E) = iVmky/l + e x 



X 2 



X 



V 



+ a/1- 



X4WR 



£ A 



WjX 



3-2 j 



(13) 

Here the square roots can be expanded using binomial series. The action variable Jcwr{E), 
defined as ^ § Cw Pcwr(x, E)dx, where Cwr is a counterclockwise rectangular contour that 
hugs the branch cut connecting Xi and x 2 , is evaluated up to order e 1 by expanding the square 
roots in Eq. ( fTBl and integrating using Cauchy's residue theorem. That yields 



J CW r{E) = — (2iri)A W2 w % . iVmfc 

Z7T 



1 + 



3 E 
16 mc 2 



— ^1 + — e 

wo V 16 



(14) 

This differs from — , the action variable for the non-relativistic harmonic oscillator, with a 
correction of order e. The angular frequency is obtained from 

1 OJcwr 



^WR 



dE 



1 m 3 , 
— ( 1 + -e) 



(15) 



which shows that the weak relativistic correction results in a fractional decrease of about |e 
in the oscillator's frequency. 

We now consider the second form of the action variable for the weakly relativis- 
tic case. The orbit equation can be solved for the coordinate to yield xcwr(p, E) = 



E-1- + 



2m 8mV 

coordinate vanishes, are defined by 



. The two physical turning momenta of this oscillator, where the 



PWRl — PWR2 



mc 



l-(l-2e) 



\ZlmE 



e 



(16) 



These are very nearly the non-relativistic turning momenta, given by — p\ = p 2 = V 2mE. 
They contain a weak relativistic correction of order e. There are two additional branch 
points, PwR3,4 °f (weak) relativistic origin given by 



PWR3 — Pwra 



We write xcwr(p, E) in the form 



V2- 



mc 



;i - 2e) 



(17) 



xcwr(p, E) = ^2 A WjP 



3-2j 



—I 



j=-oo 



'mk 



1 1 / 

yl + 7! (1 - 2£) 



P 



p 



PWRA, 
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Using this we extend xcwr(p, E) into the complex-p plane. We choose one branch cut of 
x cwr{Vi E) from piwn to P2WR along the real axis (the function is chosen positive just below 
this cut), and two other cuts, each of which joins P3WR and Pwr to p = oo along the real 
axis. The alternate definition of Jcwr(E) is 



Jcwr(E) — — — 
In 



xcwR(p,E)dp, 



(18) 



where the counterclockwise rectangular contour C' WR wraps around the branch cut connect- 
ing Piwr and P2WR- For pwR2 < \p\ < PwRi we expand xcwr m a Laurent series and 
evaluate Jqwr to obtain 

F, T 2 1 2 r 1 /t7 TI ,oo\ 2 

JcWR 



E 

LO 



1 + VT^2e 



1 _ I ( PWR2 ] 

8 \Pwra) 



(19) 



This representation of the action variable contains a series in powers of (j>wR2/PwmY 



found from 
I- dJcwR 



. To order e it is ^[1 + j^e], consistent with (JTJJ. The angular frequency is 



UWR 



dE 



1 d 

ujq de 



1 + Vl^Te 
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(20) 

which agrees with the result in Eq. ( fl5i) . Further both results agree, up to order e , with the 
expressions for the angular frequency obtained by using the fully relativistic Hamiltonian — . 



IV. HARMONIC OSCILLATOR IN QUANTUM HAMILTON- J ACOBI THEORY 
A. Quantum Action Variable 

We summarize here the formalism of Hamilton- Jacobi quantum mechanics*^, equivalent 
to other better known ones, and apply it to the special case of 1-D one particle systems 
with Hamiltonians of the form H = j- + V(x). The measurable values of the observables 
H,p and x are their eigenvalues E, p and x respectively. The equations of quantum canon- 
ical transformation are written in terms of the eigenvalues and functions of eigenvalues of 
these observables. Using the quantum characteristic function W(x, P) these transformation 
equations, equivalent to the classical ones in Eq. (jHJ), are 

dW(x,P) dW(x,P) 
P ~ dx ' A " OP • {ZL) 



The quantum Hamilton- J acobi equation, for systems with time independent Hamiltonians, 
is 

_ it fw {x ,E iP) ) + ,m(* ,e ( p))V = _ 

OX \ ox ) 

This can be either postulated or derived from Schrdinger's equation for the state ij;(x,t), 
which is written in the form e^ s ^ x,t \ where S(x,t) is Hamilton's principal function in the 
quantum context. For systems with time independent Hamiltonians, W and S are related 
through S(x,P,t) = W(x,P) — Et. Physical boundary conditions have to be imposed on 
W(x, E(P)) to complete its definition. We note that this equation resembles the classical 
Hamilton- Jacobi equation (jSJ), except for the additional term involving H, and reduces to 
it in the limit H — > 0. The dynamics described by this equation is non-relativistic and is 
equivalent in all respects to the Schrdinger formalism of quantum mechanics. We can use 
such an equation only in a perturbative sense to treat systems with weak relativistic terms 
in their Hamiltonian. 

p(x, E) = 9W ( X ^( E )) j s th e quantum analog of the classical momentum function pc{x, E) 
and, following Ref. 10, will be referred to as the quantum momentum function. Using this 
definition of p(x, E) in Eq. ( 1221) we obtain 

-ih^^-+p 2 (x,E) = 2m[E-V{x)} = p c 2 (x,E). (23) 

We note that the square of this quantum momentum function reduces to the square of 
±p c (x, E), the classical momentum function in the phase space orbit equation, in the limit 
H — ► 0. The physical boundary condition on p(x, E) is that for all x 

hmp(x,E) = +p c (x,E). (24) 
h— >o 

Eq. (T2"3"l) . along with (|24l) . defines p(x,E) which we will use in the construction of the 
quantum action variable J. Equation (123]) is the Ricatti (nonlinear) form of the Schrdinger 

iEt 

(linear) differential equation. Extending the wave function t/)^(x,t) = (f>(x, E)e~~ into 
the complex x plane it is easily shown that the nodes of <p(x, E) are also the poles of 
p(x, E). Using oscillation theorems for the linear differential equation it can be shown that 
(i) p(x, E) has simple poles of residue — ih on the real axis between the two turning points 
independent of the system's energy, and (ii) there are fixed poles whose locations and residues 
are determined by the specific nature of the potential.— 1 ^. The number of poles between 
the turning points counts the level of excitation of this quantum system. 
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Further development of this quantum formalism for periodic systems, where action-angle 
variables can be employed, requires the specification of a new momentum P, which is the 
quantum action variable J. Following the definition of the classical action variable it is 
defined as the contour integral in the complex x plane, 

J(E) = — I p(x,E)dx, (25) 
2tt Jc 

with the counterclockwise rectangular contour C tightly wrapping the real axis between the 
two physical turning points of the classical momentum function pc(x,E) = wE — V(x). 
Since C encloses only a finite number of poles of p(x, E), deforming C to encircle the poles 
and integrating using the residue theorem leads to the discretization of J: 

J= (2m){ n {-ih)} = nh. (26) 

This is in contrast to the classical action variable which assumes continuous values. At 
first glance, Eq. ff26l) seems to be the Wilson-Sommerfeld quantization condition or the 
JWKB quantum rule, but it is not. The action variable (or the Sommerfeld phase integral) 
used in those conditions is the classical action variable; the one used here is the quantum 
action variable whose definition is based on the quantum momentum function p(x, E) which 
is a solution of Eq. (I2U1) whose basis is in quantum mechanics. Finally, By deforming C 
outward and integrating we capture the energy dependence of the quantum action variable 
and obtain J(E). Thus J(E) = nh, and inverting this, we obtain the system's energy 
eigenvalues, E = E(J = nh). 

As will be shown later, the integral in Eq. (1251) can be performed without obtaining a so- 
lution of Eq. ff23l) all over the complex x plane. The quantum energy eigenvalues of a system 
that is classically periodic can thus be obtained more simply by using the discretizatized 
nature of the quantum action variable than by imposing boundary conditions on the wave 
function. It can be shown that the poles of p(x, E) on the real axis between the physical 
turning points coalesce in the classical limit, and form the branch cut of p c (x, E). The mech- 
anism of this pole coalescing is demonstrated in Ref. 10. The residue of a pole of p(x, E) 
between the turning points is proportional to h, and so is the spacing between neighboring 
poles. In the limit h — > 0, for a fixed energy of the oscillator, the real axis between the 
turning points is riddled with poles, with the quantum momentum function having opposite 
signs on the sides of the pole above and below the real axis, giving rise to the branch cut of 
p c (x,E). 
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We will construct two such equivalent forms of the action variable J for the simple 
harmonic oscillator in the quantum context and demonstrate that both forms of J are 
discrete. Since the energy is a function of the action variable, it too is discrete rather 
than continuous. It is important to recognize that (1) this quantum action variable differs 
from the classical action variable employed by Sommerfeld and others for quantization, 
(2) the quantum action variable naturally assumes values that are integral multiples of h 
and (3) the energy eigenvalues obtained by evaluating the action variable are exact. By 
making appropriate approximations of this action variable as a function of energy we obtain 
approximate energy eigenvalues. 



B. <f pdx form of quantum action variable 

The quantum action variable for the harmonic oscillator is defined through Eq. ( 1251) where 
the rectangular counterclockwise contour C closely wraps around the real axis between the 
turning points x\ and x<i- For evaluating the integral we deform C outward into the circular 
contour 7 centered at the origin, and use the Laurent series for p(x, E) in an origin centered 
annulus that includes 7. The boundary condition on p(x, E) indicates the form of the 
Laurent series that it has in this annulus. The Laurent series for the classical momentum 
function is 



Pc{x, E) = \f E — ^kx 2 = iVmkx 



We choose the form for the series for the quantum momentum function that is similar to 
that of the series for pc{x, E) and write it as 

00 

p{x,E) = Y,bjX 3 ~ 23 . (28) 
3=1 

Substituting this in Eq. (1231) and equating coefficients of like powers of x we get b\ = ±Wmk. 
We choose the + sign here since ai = +i\ / mk; the boundary condition on p{x, E) requires 
bj — > dj in the classical limit. The next coefficient is 62 = — ^a/t^ + *f • Using the residue 

theorem we get 

J ( E) = 2,i(±) h = E^- \ (29) 

This, along with the discretization condition in Eq. (1261) . yields E = (n + ~)H, the energy 
eigenvalues of the harmonic oscillator. 
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C. Equivalent quantum Hamilton- Jacobi equation 



We present here a new version of the quantum Hamilton- Jacobi equation that is equivalent 
to the one in Eq. (|22p using an alternate quantum canonical transformation scheme. While 
the dynamical equation can be postulated for the quantum generating function S(p, X, t) 
that generates the transformation 



x 



dS(p,X,t) 



p 



dS(p,X,t) 



dp dX 

and show its equivalence to the Schrdinger equation, we choose the reverse route and derive 
it from the latter. Starting with the Schrdinger equation 



> 



ih 



d\<f> > 



(30) 



and using the momentum representation where x — > ih-^-, p — > p and \<f> >— > (j){p, E, t), we 
get 



p 



— + v { ih— 



2m 



d 



dp 



iH M- 



(31) 



Introducing the alternate quantum Hamilton's principle function S(p,X,t) through 
4>(p, E(X),t) = e i s( - p ' x '^ and defining G from 



d 

\/(^— ) (j)(p, E,t) = G[h, S p , Spp, . . . 



eh 



S(p,X,t) 



we obtain 



^— + G(h, Sp, Spp, . . .) 



-St. 



(32) 



(33) 



Here suffixes denote partial differentiation. We define the alternate quantum characteristic 
function W(p, X) through the relation S{p,X,t) = W(p, X) — Et and use it in Eq. (|33|) to 

get 

E. (34) 



P 



j^ + G(h,W p ,W P p,... 

This is the quantum Hamilton- Jacobi equation for W(p,X) which generates a canonical 
transformation from (x, p) to (X, P). Use of the quantum canonical transformation equation, 



x = — 



aw 

dp ' 



in Eq. (G£Q) leads to 



— + G[h, x{p, E), . 

2m \ op 



E. 



(35) 
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This is the equivalent of Eq. (1231) under this alternate canonical transformation scheme. 
The physical boundary condition on x(p,E) is that in the limit H — > 0, it should reduce to 
the classical coordinate function, xc{p,E). Specializing to the case of the simple harmonic 
oscillator we get 

(36) 



q = i dxjp, E) _ }_ x 2, p E ^ 

h dp h 2 ' 



Substituting this in Eq. fl35|) we obtain the quantum equivalent of the classical orbit equation: 



Bx 2 2 
in— + x = - 
op k 



E 



p* 
2m 



(37) 



In the limit h — > this equation becomes, for real values of x and p, the classical orbit 
equation that constrains the harmonic oscillator to an elliptical path in phase space. 

We choose that quantum canonical transformation which makes the quantum action 
variable J the new coordinate. It is defined as the contour integral 



J = —^-<j>^x(p,E)dp 



(38) 



in the complex-p plane, with the clockwise rectangular contour C enclosing the real axis 
between the two turning momenta p\ and p2- A comparison of Eq. ( 1371) with Eq. ( 123]) shows 
that, as in the case of the quantum momentum function p(x,E), the quantum coordinate 
function x(p, E) has a finite number of poles on the real axis between p\ and P2, each of 
residue +ih. In the classical limit these poles coalesce to form the branch cut of xc{p,E). 
The steps for the determination of energy eigenvalues here parallels that described in the 
previous section. We evaluate J in two ways, one by outward deformation of C into the 
origin centered circular contour 7' and integrating to obtain J = J{E), and, two, through 
inward deformation, encircling the poles and integrating, to get J = nh. 

For \p\ > p 2 , the form of the Laurent series for the classical momentum function xc{p, E) 



is 



x c (p, E) 



E — 



2m 



'mk 



P 



3-2j 



j=-oo 



We use a similar series form for the quantum coordinate function x(p,E) valid for this region 
of the complex plane: 

00 

x(p,E)=J2b' 3 p 3 - 2j . (39) 



14 



Substituting this in Eq. (I35I) . and imposing the physical boundary condition on x(p, E), we 
obtain 

b l = -7=, &2 = -* I o 

ymk \^ 
Expanding the contour C outward to 7' and evaluating J we get 

J = - -L (27ri) 6' 2 = - - \ (41) 

Z7T CJo ^ 

Shrinking the contour inward to encircle the poles on the real axis between p\ and P2 and 
integrating yields J = nh. So, ~ — f = and inverting this relation we obtain the energy 
eigenvalues of the simple harmonic oscillator, E — (n + |)^o- 

This alternate form of the quantum action variable has the same physical content as the 
previous form. While the first form is easier to evaluate for systems whose potential energy 
functions are not necessarily quadratic in x, the second form is better suited for the weakly 
relativistic oscillator, as will be shown in the next section. 




V. QUANTUM MECHANICS OF WEAKLY RELATIVISTIC OSCILLATOR 



A. Integral <f pdx form of action variable 



The Hamiltonian for the weakly relativistic oscillator is j- — Sl ^a c i + \kx 2 . Using this 
in the quantum Hamilton- Jacobi equation for W(x, P) (we omit the subscript WR for W), 
which is obtained from Schrdinger's equation following the steps outlined in Section IV. A, 
we get 

2 d 2 W 



E 



-kx 2 
2 



1 fdW 



2m \ dx 



ih d 2 W 
2m dx 2 



1 

3m 3 c 2 



dW 
dx 



6ih 



dWd 3 W o fd 2 W 
h 2 < 4— — + 3 



dx dx 3 



dx" 2 




(42) 



Using the quantum canonical transformation equation p = dW/ dx in Eq. (T42|) results in 



E — 
h 



P 2 I 1 U X 2 P 4 

2m 2 8m 3 c 2 



ih dp 
2m dx 2m 3 c 2 



3 2 dp I 3 / dp 



+ p 



d 2 p 1 ih 2 d 3 p 



dx 2 



4 dx 3 



(43) 
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Eq. (|4"3l) shows the presence of terms proportional to h or its higher powers, which are absent 
in the corresponding classical orbit equation. Further, in the non-relativistic limit it reduces 
to Eq. (j23j). 

Using the solution Pwr of Eq. (T4"51) satisfying the physical boundary condition (i.e., has 
the correct classical limit) we define the quantum action variable as 



WR 



2tt 



— <p pwR(x,E)dx, 



(44) 



using the same contour Cwr as in the classical case. Following the form of Pcwr(%, E) in 
Eq. ( |T3l) the Laurent series for Pwr{%, E) in the annulus x<i < \x\ < x^wr, correct to order 
e 1 , is of the form 

2' 



p WR (x,E) 




(45) 



The non-relativistic quantum momentum function in Eq. ([28]) has been modified here by 
the addition of a relativistic term of order e. We have chosen the form of this series such 
that (i) the relativistic modification is of the first order in e and in the limit e — ► we 
get Pwr(x, E) — > p(x, E), and (ii) in the limit h — » we obtain Pwr{x, E) — > pcwr{%, E). 
The coefficients bj are known from the non-relativistic quantum case previously considered 
(see Section IV. B). Substituting this series form of the solution in Eq. (143j) we obtain the 
coefficients B and Bi\ 



Bn 



1, B 1 = l + 



~AE 



(46) 



The coefficient of x~ x in Pwr(x, E) is imuo x\ (biB — b Bi). Evaluating J by deforming 
the contour Cwr outward to the origin centered circular contour ^wr previously considered, 
we get 



WR 



E 

UJ Q 



1 




imuo x\ (bxBo - 6 ^i) 

2~ 



(47) 



A comparison with Eq. ( 1291) shows the relativistic corrections present in this quantum 
action variable up to order e. Solving Eq. ( 1471) for E, we get 

2 



E = E(J WR = nh) 



n 




n + 




huJc] 



(4£ 
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This gives energy eigenvalues with a first order relativistic correction that is proportional 
to the dimensionless energy parameter which measures the level separation in the 
non-relativistic case. The energy levels in the weak relativistic case are lower than in the 
non-relativistic case, and the lowering is predominantly quadratic in the quantum number 
n. The seperation between energy levels is proportional to n. 



B. — <f xdp form of action variable 

We now apply the formalism developed in Section IV. C to construct the quantum action 
variable in its alternate form. Using the Hamiltonian for the weakly relativistic oscillator 
in the quantum Hamilton- Jacobi equation for the alternate characteristic function W(p,X) 
we get 



2 



B 2 W ldW\ 2 

— in- 



2m 8m 3 c 2 



(49) 



dp 2 \ dp J k 

In comparison to the quantum Hamilton- Jacobi equation ( 1421) . this alternate form is simpler 
as it has a single quantum term, and is to be preferred in the study of this oscillator. 
Introducing the quantum canonical transformation equation x = —dW /dp in Eq. ( 1491) . we 
obtain 



dx 2 2 



ih— + x 
op k 



2m 8m 3 c 2 



(50) 



We impose the physical boundary condition that in the limit h — > 0, the quantum coordinate 
function x(p,E) should reduce to the classical coordinate function xc(p,E) for all p. This 
equation has the same structure as its non-relativistic counterpart, Eq. (1371) . 

Following the classical case, we define the quantum action variable as Jwr — 
— 2^ic' x wr(Pi E)dp, where xwr is the solution of Eq. (150]) that satisfies the physical 
boundary condition. The construction of this form of the quantum action variable and the 
derivation of energy eigenvalues in this case is identical in all relevant details to the problem 
of the the non-relativistic anharmonic oscillator with a quartic potential energy term Sx , 
and with J in the § p(x, E)dx form, shown in the Appendix. A comparison with the latter 
problem indicates that we can find the action variable for the weakly relativistic oscillator 
by making the replacements x — > —p,p — > —x, -p — > in Eq. (1571) and get 

T , E h 3h f 4E 2 } fhw \ 
JwR{E) = - - 2 + 64 ( 1+ (^f) Wj' (51) 
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correct to the first order in Solving for E in this order using Jwr(E) = nh, we get 

2 



E 



n + 




n + 



1 



+ 4 



(52) 
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The above energy eigenvalues obtained by approximating the quantum coordinate func- 
tion, while not identical to those in Eq. f j48l) . display the same features that were observed 
earlier. Both forms of the quantum action variable indicate that the lowering of energy levels 
due to the weak relativistic correction is approximately proportional to n 2 . 



VI. OSCILLATOR FREQUENCIES AND ENERGY LEVELS FROM DIFFERENT 
SCHEMES: A COMPARISON 



Table 1 summarizes the expressions for the classical action variable of the relativistic 
oscillator under four different schemes, two of which we have considered here. The first two 
rows show two equivalent series representations for Jc for the fully relativistic oscillator with 
no approximation^, and the next two rows the corresponding representations obtained by 
using the weak relativistic approximation for the kinetic energy. To the first order in e they 
all yield the same expression. Relativistic dynamics lowers the oscillator's frequency due 
to time dilation, and the fraction by which the frequency is lowered is |e for this weakly 
relativistic oscillator. The quantum action variable for the weakly relativistic oscillator under 
the two alternate action variable schemes we considered is shown in Table 2. The energy 
eigenvalues obtained from these schemes and those from two other approximation schemes 
are also shown in this table. The semiclassical JWKB approximation involves discretization 
of the classical action variable Jcwr(E) in Eq. (fTi]) for the weakly relativistic oscillator: 



J, 



CWR 



E 

LO 



16 



(n + -) h. 



In the Rayleigh-Schrdinger perturbation scheme the shifts in energy eigenvalues from their 



p 



8m 3 c 2 



unperturbed values, to the first order in e, are AE n =< <p n \H e \(j) n >, treating H e - 
as the perturbation term in the Hamiltonian whose zeroth order form is Ho = $~ + ^kx 2 . 
In all four schemes the calculated weak relativistic shift in the energy eigenvalues from thier 
non-relativistic values is approximately — ^ (^t) n2 - We find that the energy level spacing 
is 



E n+1 - E n w 



3 fUjjQ 

8 mc z 



3 

1 - - e 



ftw 



WR- 
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The angular frequency ujwr, unlike in the non-relativistic case, is energy dependent. Thus 
the spacing between the quantum energy levels in a range which conform to this approxima- 
tion scheme is proportional to the classical oscillator's weak relativistic angular frequency for 
that energy range. Viewed differently, the weak relativistic correction produces a fractional 
shift of |e in the angular frequency of the classical oscillator. There is an identical fractional 
shift in the energy level separation of the corresponding quantum oscillator. 

VII. QUANTUM - CLASSICAL AND RELATIVISTIC - NON-RELATIVISTIC 
CORRESPONDENCE 

The physically correct theory of matter should be both quantum in nature and meet the 
principle of relativity. The Klein-Gordon and Dirac equations represent physical theories 
that meet these two physical requirements. The description of the weakly relativistic oscil- 
lator considered here meets the first principle. While it does not meet the second principle 
it does incorporate a weak relativistic dynamical correction. It is a useful treatment of an 
"energetic" oscillator as it begins to approach the relativistic regime. We notice the op- 
eration of two correspondence principles in the dynamics of this physical system. One is 
the reduction of the relativistic model to its appropriate, well established non-relativistic 
form. The other is the similar reduction of the quantum model to its corresponding classical 
form. Both the weak relativistic quantum Hamilton- Jacobi equations ( )42l) and ( H9j) reduce 
to their weak relativistic classical counterparts in the limit h — > 0, and these equations, 
in turn, assume their classical non-relativistic form for e << 1. Secondly, the poles of the 
quantum momentum function Pwr{ x , E) in the complex-x plane, and those of the quantum 
coordinate function xwr{p, E) in the complex-p plane coalesce to form the branch cuts of 
classical functions Pcwr( x , E) and xcwr(Pi E) respectively. Thirdly, the expressions for 
Jwr(E) in the quantum case reduce to those of Jcwr(E) in the classical non-relativistic 
limit. Fourthly, the weakly relativistic quantum action variable Jwr{E) reduces to the non- 
relativistic quantum action variable J(E) for e << 1. This formalism of quantum mechanics, 
applied to the harmonic oscillator, thus shows in a unified framework the critical role played 
by the two "small" parameters, h and e, in the emergence of non-relativistic classical theory 
from non-relativistic quantum theory and a non-relativistic quantum model from a weakly 
relativistic quantum model. 
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VIII. CONCLUSION 



The use of action-angle variables in the study of classical periodic systems has a parallel 
in quantum mechanics. Using the quantum action variable the exact energy eigenvalues of a 
bound system that is classically periodic can be obtained without a detailed solution of the 
dynamical equations. We have demonstrated the construction of two forms of the quantum 
action variable, and applied it to the harmonic oscillator to obtain its energy levels. We have 
extended the use of this formalism to a weakly relativistic harmonic oscillator. The classical 
frequency of such an oscillator is lowered from its non-relativistic value by the fraction |e. 
There is an identical fractional shift in the energy level separation of the corresponding 
quantum harmonic oscillator due to the leading order relativistic correction. While the 
problem must be properly addressed by a relativistic quantum theory we have shown how 
this system may be studied in a unified manner within a non-relativistic quantum Hamilton- 
Jacobi theory. 

IX. APPENDIX: ANHARMONIC OSCILLATOR WITH QUARTIC TERM 

We consider an approximation scheme for obtaining the energy eigenvalues for a quartic 
anharmonic oscillator described by the potential V(x) = \kx 2 + <5x 4 , with 5 > 0, using the 
quantum action variable. The 5 < case can be treated in a similar manner. The term 5x 4 
is considered "small" for \x\ < x 2 ^ in comparison to the dominant term |/cx 2 . Dynamical 
variables specific to this oscillator have the suffix AHO. We denote the physical turning 
points for the 5 = (or the simple harmonic oscillator) case by —xf^ = x^ = #i 
and Xi are the physical turning points in the presence of the quartic potential term, and are 
approximately of mag nitude xf (l - f|<5). For 5 > they are closer to each other than 
those in the harmonic oscillator case. We seek the energy eigenvalues of this anharmonic 
oscillator correct to order S 1 . The classical momentum function pcaho{x, E) has two other 
(unphysical) turning points, approximately given by — x 3 = x 4 ~ i \J~^ (l + jjr)- We write 
the classical momentum function in the form 



2 




(0) 



) 



2 



2 



Pcaho(x,E) = iVmk x 1 



1 + 




) 



2 



(53) 
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To order 5 1 this has, for x 2 < \ x \ < l x 4|) the series representation 



Pcaho(x,E) 



mk x 




, 5 

1 + T x 
k 



oo 

E 

2=0 



(0)' 



2/ ' 



(54) 



Here Cj are the coefficients in the binomial expansion of y/1 — u in powers of u. We use a 
similar structure for the Laurent series, for the anharmonic oscillator's quantum momentum 
function paho{x, E), good for x 2 < \x\ < \x^\ and correct to order 5 1 : 



Paho(x,E) 



mk x 




, 5 

1 + T x 
k 



1=0 



x 



(0)' 



21 



X 



(55) 



The coefficients bj are known from our solution of the quantum harmonic oscillator problem 
(Section IV. B). We need to solve for the coefficients D\ using Eq. fl23|) . The quantum action 
variable is defined as 

J AHO {E) = ^-(p p AHO {x,E)dx, (56) 



Cai 



where the counterclockwise rectangular contour Caho wraps around the real axis between 
the turning points x\ and X2 (or xf* and x 2 for the 5 < case). The coefficient of x~ l in 
the series in Eq. ( 1551) is imoj^x^ &i + f xf^ (b 2 D + boD 2 + biD\) . 

We define the dimensionless parameter A = which is "small" for large values of the 
oscillator's energy. Introducing the above series for paho{x, E) in Eq. ( 1551) and solving for 
the D-coefficients, we get D = 1, D x = 1 + A, £> 2 = 1 - |A + 2A 2 . 

Evaluating Jaho by the outward deformation of the contour Caho m to the circular 
counterclockwise contour ■Jaho we obtain 

11 38H /. 1 
" + 4A 



(57) 



2 ' 4A 32m 2 ^ V' ^ 2 
Solving for E in Jaho{E) = nh (obtained from inwardly deforming Caho to encircle the 

poles and integrating) we get the energy eigenvalues of the quartic anharmonic oscillator, 

to order 5 1 : 



E = fkv, 



o 



1 3/5 



TlJjJc 



hujQ < n + 



1 



3/5 
8 \T? 



huJr 



(58) 



For 5 > 0, the quartic perturbation shifts the oscillator's energy levels higher (and lower for 
the 5 < case), with the spacing between the energy levels increasing linearly with n. 
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Tables 



TABLE I: Action variable for a classical relativistic harmonic oscillator from 4 schemes. 



Nature of oscillator Form of Jq Expression for Jc{E) 


Fully relativistic case J- § pdx 


+ D i-l(^)-^(^) 2 --- 


Fully relativistic case xdp 


£aA+I [i-^ + ^ + t^ 3 ...] 


Weak relativistic case ^- <f pdx 


^ (1 + ^e), to order e 


Weak relativistic case xdp 


f (1 + to order e 



TABLE II: Correction to zeroth order energy eigenvalues, E n $ = (n+ ^) /io>o, of weakly relativistic 
harmonic oscillator. 



Calculation scheme Calculated term Correction to K 



n,0 



Evaluate quantum action J WR = JL[1 + C {£ + -^M) [{n + |) 2 - (|) 2 ] 



variable in ^ § pdx form jq 


C frw 
\ E J 






Evaluate quantum action Jwr = ^ — § 




j- + i-!(^)} 2 + 4 


mc 2 



variable in f sdp form + j {l + ^} ffi 



JWKB approximation J CW R = J_[l + ^e] ~^^o(n + ±) 2 ^ 

Rayleigh-Schrdinger A£„ =< 0„| - ^|0 n > -^o[(n + ±) 2 + |] ^ 

perturbation theory 
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